Stability of cnoidal waves in the parametrically driven nonlinear Schrodinger equation 
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The parametrically driven, damped nonlinear Schrodinger equation has two cn- and two dn-wave 
solutions. We show that one pair of the cn and dn solutions is unstable for any combination of the 
driver's strength, dissipation coefficient and spatial period of the wave; this instability is against 
periodic perturbations. The second dn-wave solution is shown to be unstable against antiperiodic 
perturbations — in a certain region of the parameter space. We also consider quasiperiodic per- 
turbations with long modulation wavelength, in the limit where the driving strength is only weakly 
exceeding the damping coefficient. 
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I. INTRODUCTION 

Periodic waves arise in media of different physical na- 
ture. They find direct applications in plasma physics |l|, 
2], nonlinear optics solid state physics [1,13], an( ^ 

physics of Bose- Einstein condensates Q. The nonlinear 
periodic waves (or cnoidal waves) interpolate between 
plane waves and solitons. 

One of the most important problems associated with 
periodic waves is their stability in various nonlinear me- 
dia. To date, stable periodic patterns were discovered in 
physical settings modelled by self-defocusing nonlinear- 
ities and in Bose-Einstein condensates confined 

by periodic potentials [l4|, EB| • Families of stable cnoidal 
waves supported by focusing nonlinearities were revealed 
in quadratic media [H, [ljj ■ All these results pertain to 
conservative media where self-supported structures exist 
due to the balance between dispersion and nonlincarity. 
In the presence of dissipation, an additional balance be- 
tween gain and loss is needed. In the dissipative case the 
parameters of the solution (e.g., amplitude and phase) 
are fixed by the parameters of the governing equation, 
whereas in conservative settings solutions form continu- 
ous families. So far, stable stationary dissipative periodic 
waves were mostly observed in optical cavities. 

In this paper we study stability properties of periodic 
wave solutions of the parametrically driven, damped non- 
linear Schrodinger equation: 

i4> t + ^xx-i> + 2i)\i)\ 2 = hvb* -i^, h,j>o. (l.i) 

Here 7 is the damping coefficient and h the amplitude 
of the parametric driver. Equation (|1.1[) is an archety- 
pal equation for small and slowly-varying amplitudes of 
waves and patterns in spatially-distributed parametri- 
cally driven systems. It arises in a wide variety of physical 
contexts including instabilities in plasma |22|. |23l . am- 
plitude generation in Josephson junctions 124 , 125I 1. and 
signal amplification effects in fiber optics 
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At the same time, the knowledge about periodic non- 
linear waves of PDNLS is rather scarce. Umeki [2^] ex- 
amined numerically stability properties of the families of 
cn- and dn- waves of PDNLS considered in the context of 
water waves in a vertically forced long container {h < 0). 
For a fixed value of the spatial period of the wave, he 
obtained the corresponding stability diagram and stud- 
ied the temporal evolution of the perturbation. Umeki's 
numerical considerations were restricted to perturbations 
of the initial wave profile having the same period as the 
cnoidal wave itself. 

In this paper, we examine the linear stability of the cn 
and dn type solutions of Eq. (|1.1[) with respect to two 
broad classes of perturbations. The first class consists of 
perturbations which are periodic or antiperiodic with the 
period of the wave. These perturbations have a discrete 
spectrum and are relatively easy to analyse. On the other 
hand, the periodic or antiperiodic perturbations turn out 
to be sufficient to detect the instability of the underlying 
wave in a wide range of parameter values. 

In particular, we prove that out of two cn- and two dn- 
wave solutions, one pair of cn- and dn-waves is unstable 
against periodic perturbations — for any combination of 
the driver's strength, dissipation coefficient and spatial 
period of the wave. The other dn wave is shown to be 
unstable against antiperiodic perturbations — in a cer- 
tain region of the parameter space. 

The second class includes quasiperiodic perturbations 
which have the form of symmetry eigenvectors modulated 
by a long wavelength. Here, the small wavenumber of the 
modulation serves as a perturbation expansion parame- 
ter. The perturbation theory can be developed near the 
lower boundary of the cnoidal-wave existence domain, 
that is, for driving strengths only slightly exceeding the 
damping coefficient. 

Using this perturbation approach, we reproduce sev- 
eral results obtained by other means. First, we confirm 
the instability of the cn- and dn-wave solutions of the 
unperturbed (7 = h = 0) NLS equation, a result already 
available in literature. Second, we corroborate our own 
conclusions on the instability of the two damped-driven 
dn waves in the h ss 7 case. More importantly, the long- 
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wavelength-modulation treatment provides useful infor- 
mation on the structure of the spectrum of the cn wave 
whose stability cannot be classified by restricting to the 
periodic perturbations. 

We complement the perturbation study of the cn wave 
with a numerical analysis of its linearised spectrum in 
the h ss 7 case. 

II. PERIODIC SOLUTIONS AND THEIR 
LINEARIZATIONS 

A. Two pairs of cnoidal waves 

Eq. (|l.ip has two pairs of nonequivalent periodic so- 
lutions. One pair is expressible in terms of the Jacobi 
cosine function: 

V? = A ± q cn (A ± x,k)e-^, (2.1) 

where 

1/V2 < k < 1. We will occasionally be referring to these 
solutions as cn + and cn~, respectively. The other two 
solutions (referred to as dn + and dn - in what follows) 
invoke the Jacobi dn function instead: 

*i n = A±q dn (A±x,k)e- ie ±, (2.3) 

where 

*<° {x ' k) = v^Hvrh- h )- (2 - 4) 

< k < 1. The amplitudes A± are given, in both cases, 
by 

A 2 ± = l± Vh 2 -J 2 , (2.5) 

and the phases by 

1 /7\ n 

9+ = -arcsin — , = 6>. 

+ 2 KhJ 2 + 

Eq. (|2.5[) carries the entire information on the domain 
of existence of the four periodic solutions on the (/i/y)- 
parameter plane. For the given 7, the cnoidal waves "J™ 
and 'J'^" exist for all h > 7 whereas the domain of the 
solutions ^1" and is bounded on both sides: 7 < h < 

y/l + 7 2 . The two cn solutions have the spatial periods 
L C n/A±, where 

Lcn = L cn (k) = 2K{k)^2k 2 - 1, (2.6) 

and the dn solutions are periodic with the periods 
Ldn/A±, where 

idn = L dn (k) = 2K(k)^2-k 2 . (2.7) 



In Eqs. (|2.6l ) and (|2.7p . K is the complete elliptic integral 
of the first kind. As k varies from l/\/2 to 1, L cn grows, 
monotonically, from to infinity. Less obvious fact is that 
L^n is also a monotonically growing function - growing 
from \Z2tt to infinity as k varies from to 1. (See the 
Appendix.) Therefore, for the given h and 7 the period 
of the cnoidal wave can be used as its third parameter, 
in lieu of the elliptic modulus k. 

As k — > 1, the periodic solutions ^P™ and ^^j" approach 
the tb+ soliton, while the $™ and <J/! n tend to the soliton 

i>- El- 

B. Linearised problem 

To examine the stability of cnoidal waves (|2.ip - (|2.3p . 
we write ip(x,t) — 1 &±(x) + 6ip(x,t), where ^±(x) is the 
stationary solution in question, and linearize in 5tp. Let- 
ting 

Sip(x, t) = e- z9± [U(x, t) + iV(x, t)} 

yields 

U T = C Q V, -V t -2*/V = CiU, (2.8) 
where the operators £0,1 are given by 
d 2 

£ = --^ + lT£-2q 2 (2.9a) 

and 

d 2 

A = -^2+l-6g 2 . (2.9b) 

In (|2.8p . T — A\t is the scaled time variable, and in 
(I2.9ap - (l2.9bp . X — A±x is the scaled spatial coordinate. 
In (|2.8p ~ (|2.9p we have scaled the damping coefficient and 
introduced a parameter £ which measures the driving- 
damping difference: 

7 = 7M±, £^2^h?^/A 2 ± . (2.10) 

The notation q stands for q cn or q^n [Eq. (|2.2p or 
Eq. (|2.4p ] depending on whether we linearize about a cn 
or dn solution. The top sign in front of £ in Eq. (|2.9ap 
corresponds to the and \?^ n solutions; the bottom 
sign selects the ^P™ and \l/l n cnoidal waves. The cnoidal 
wave is deemed unstable provided Eqs. (|2.8p have solu- 
tions growing faster than exp(7T) in time. 

We do not impose any periodicity requirements on U 
and V. All we require is that U(X, T) and V(X, T) be 
bounded on the whole line — 00 < X < +00. 

We call A a point of spectrum of an operator C if the 
equation Cy = Ay has a solution y(X) bounded for all 
X, 00 < X < 00. The spectrum of the operators (|2.9ap 
and (|2.9bp can be found exactly. Consider first the cn 
solutions, that is, assume that q(X) in Eqs. (I2.9a|) ()2.9bp 
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is given by Eq. ([2"T2|) . Then, denning £ = X/^2k 2 - 1, 
we obtain: 



and 



where 



. _ £< 2 )-l-2fc 2 



2fc 2 - 1 



(2.11) 



(2.12) 



Ml) _ 



+ l)fc 2 sn 2 (£,fc) 



is the ^-gap Lame operator [30] ■ Therefore the spectra of 
Co and C\ result from the spectra of £W and £^ 2 ^ (given 
in |30|) by shift and scaling. 

Thus, the spectrum of Co consists of a finite band 



Ae 



l 



2fc 2 - 1 
and a semi-infinite band, 
fc 2 

A e 



t£, t£ 



t£, oo) . 



(2.13) 



_2fc 2 - 1 

The spectrum of C\ comprises two finite bands 

2fc 2 



-1 - 



2VT 



fc 4 



2k 2 - 1 ' 2fc 2 - 1 
and a semi-infinite band 



0, 



2(1 



2fc 2 - 1 



-1 



Vl - k 2 + fc 4 
2fc 2 - 1 : 



In the case of the dn solutions (|2.3|) the spectrum struc- 
tures are similar. Namely, assuming that q(x) is given by 
Eq. (|2.4[) and defining £ = X/ V2 — k 2 , the operators Co 
and £i can be expressed as 



and 



C 



& 1 



2-fc 2 



£( 2 



2-fc 2 



2-fc 2 



(2.14) 



(2.15) 



respectively. Accordingly, the spectrum of Co comprises 
two bands, 



i 



t£ 



i 



T^, oo , 



while the spectrum of C\ consists of three: 



-1 - 



2V1 - fc 2 + fc 4 
2-fc 2 



'2-fc 2 
-1 



2(1 -fc 2 
2-fc 2 



vr 



fc 4 



(2.16) 



-, 
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(2.17) 



In Eqs. (|2.13p - (l2.17l) the top sign in front of £ pertains 
to the ^-i- solutions and the bottom sign to 1 J r _. 

The spectrum of the ip+ and ^_ solitons (2ljj is re- 
covered be sending fc — > 1. In this case, each finite 
band collapses into a discrete eigenvalue. As a result, 
the spectrum of Cq consists of a single discrete eigen- 
value Ao = ^f£ (and a continuum of values A > 1 =p £ ) 
whereas the spectrum of C\ includes two discrete eigen- 
values, Ao = —3 and Ai = ( and a continuum A > 1). 



Stability eigenvalues and symplectic eigenvalue 
problem 



We assume that the linear system (|2.8|) has separable 
solutions of the form 

U(X, T) = Re [e^uiX)] ; V(X, T) = Re [e 7lT i(X)] , 

where u, v and rj are complex. Here r\ and u, v are 
eigenvalues and eigenfunctions in the eigenvalue problem 



C v = r/u, C\u = —(rj + 27)? 



(2.18) 



We will be referring to this problem as the "linearised 
eigenvalue problem" , whereas r\ will be called "stability 
eigenvalues" below. 

Making a substitution (u, v) — > (u,v) where Xv — rjv, 
u = u, and 



A 2 =77(77 + 27), 



we transform (|2.18p to 



where 



and 



H = 



d 



Co 



(2.19) 
(2.20) 

(2.21) 



One advantage of this formulation is that eigenvalues A 
depend on h and 7 only in combination £ (that is, only 
as h 2 — 7 2 ); thus a two-parameter eigenvalue problem is 
reduced to a one-parameter problem. Another merit is 
that the operator J~ lr H is symplectic, that is, generates a 
hamiltonian flow. The spectrum of symplectic operators 
consists of pairs of opposite pure-imaginary values, real 
pairs and complex quadruplets. If A is a real or pure 
imaginary point of spectrum, then —A is another one; if 
a complex A is in the spectrum, then so are —A, A*, and 
—A* [3lJ. We will be referring to the eigenvalue problem 
(|2.20p as the "symplectic eigenvalue problem", and A's as 
the "symplectic eigenvalues" . Having found a symplectic 
eigenvalue A, we can readily recover the corresponding 
growth rate Re 77 from (|2.19[) : 



Re 77 



-7 + Re V7 2 + A 2 . 



(2.22) 
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(Here we have kept the largest of the two growth rates.) 

Before proceeding to the analysis of the symplectic 
spectrum, three remarks are in order. First of all, 
we note that if the symplectic eigenvalues ±A are real, 
the corresponding stability eigenvalues are real as well: 
rj = —7 ± \/j 2 + A 2 . Thus the occurrence of a real sym- 
plectic eigenvalue immediately implies instability of the 
underlying cnoidal wave. This instability has monotonic 
growth. 

Second, since the potentials in the operators £0 arL d 
£1 are even, it is sufficient to consider only even and odd 
eigenfunctions y{X). 

Finally, we note that the potentials of the operators 
£0,1 are periodic with the period L cn in the case of the 
cn ± solutions and in the case of the dn ± cnoidal 
waves. Therefore, the subspace of periodic functions 
with period L cn respectively L^ n is invariant under the 
action of the operators £0,1 in the case of cn* respec- 
tively dn* solutions. This implies that the eigenvalue 
problem (|2.20[) is well posed on the subspace of peri- 
odic functions with the corresponding period. The sub- 
space of antipcriodic functions, that is, functions satisfy- 
ing y(X + L C n,dn) = —y(X), is also invariant. Accord- 
ingly, the eigenvalue problem (|2.20p is well posed on the 
subspace of antiperiodic functions. 



with the mixed boundary conditions 



L 



yx(o) = o-, y[^)=o. 



(3-3) 



One eigenvalue of £0 is Ai = =f£; it is associated with 
the eigenfunction 



2/1 ( 



(X) = cn {X/^2k 2 - 1, k 



(3.4) 



The eigenfunction does not have zeros inside the inter- 
val (0,L/2); therefore Ai is the lowest eigenvalue of the 
regular Sturm-Liouville problem (I3.2j) + (|3.3p . If we are 
considering the linearization about the solution cn~ , the 
eigenvalue Ai = +£ is strictly positive. Hence the oper- 
ator £0, defined by the differential expression (|2.9a[) and 
boundary conditions (I3.3[) . is positive definite. On the 
subspace of functions satisfying (|3.ip . the system (12.20)) 
can be written in the form 



C\u = — A 2 £q u, 



(3.5) 



where £1 is symmetric and £q 1 symmetric and positive 
definite. The smallest eigenvalue (— A 2 )o of the gener- 
alised eigenvalue problem (|3.5[) is given by the minimum 
of the Rayleigh quotient 



III. INSTABILITY TO PERIODIC AND 
ANTIPERIODIC PERTURBATIONS 

To establish instability of a solution it is sufficient to 
demonstrate its instability to a particular class of per- 
turbations. In this section we will show that the cnoidal 
waves cn~ and dn~ are unstable w.r.t. periodic pertur- 
bations, for any combination of the parameters h, 7 and 
k. We will also show that the wave dn + is unstable to 
antiperiodic perturbations — in some region of the pa- 
rameter space. 



(-A 2 ) = min 



(u\£i\u) 



« (u\C^\u) 



(3.6) 



where the minimum is evaluated over all functions u(X) 
satisfying the boundary conditions (|3.1[) . Here the scalar 
product is defined by 



L/2 



u(X)v(X)dX. 



Turning to the operator £1, we note that it has a nega- 
tive eigenvalue ji = —2k 2 / (2k 2 — 1) with the eigenfunction 



A. Instability of cn 

In the case of the cnoidal wave cn~ , we will show 
that the symplectic eigenvalue problem f|2 . 20[) has real 
eigenvalues A associated with the eigenfunctions satisfy- 
ing boundary conditions of the third kind, 

ux (0) =«(§)= °' vx (0) = u (I) = °- 

Here L = L cn is the half-period of the cnoidal wave, 
defined by Eq. (|2.6p . Due to the boundary condition at 
the origin, the eigenfunction y(X) has to be even. 

Let £0 be defined by Eq. (|2.9al) with q = q cn and con- 
sider an eigenvalue problem 

C y = Ay (3.2) 



zi(X) = ca( , X ,fc)dn( - X =,k] (3.7) 
U ' W2fc 2 -1 J W2fc 2 -1 J 

which satisfies the boundary condition (|3.3p . Therefore, 
the quadratic form (u|£i|u) attains negative values on 
the space of functions with the boundary condition (|3.3p . 
Eq. (|3.6p implies then that the smallest eigenvalue — A 2 of 
the eigenvalue problem (|3.5p is negative and hence the 
vector eigenvalue problem (12.20)) has a real eigenvalue A. 
By (|2.22p we conclude, eventually, that there are pertur- 
bations with positive growth rates Re rj. 

Since the eigenfunction y(X) is even, y(L/2) = im- 
plies y(—L/2) = 0. On the other hand, the derivative 
yx(X) is odd; thus we have 




(3.8) 



5 



Considering Eq. (|2.20[) as a system of four first-order 
equations with L-periodic coefficients, Eq. (|3.8l) implies 
that its real solution (it, ux, v, Vx) is antiperiodic on the 
interval (— L/2,L/2) — that is, has the Floquet multi- 
plier — 1. By the Floquet theorem, this solution (and 
hence the eigenfunction y) is periodic with the period 
2L. 

We conclude that the cnoidal wave cn~ is unstable 
against periodic perturbations of period 2L CD (the period 
of the wave), for any combination of h, 7 and k. 



B. Instability of dn 

The case of the wave dn~ is different from the situ- 
ation in the previous section in that the operators Cq 
and C\ have potentials periodic with the period of the 
cnoidal wave (rather than with the half period). Nev- 
ertheless, the above proof remains in place, with some 
modifications. 

In this case, the spectrum of the operator Cq lies on the 
positive semiaxis. Therefore, the operator Cq is positive 
definite for any boundary conditions imposed on y(X) in 
(|3.2p . in particular for the boundary conditions 



(3.9) 



Therefore, on the subspace of functions satisfying the 
Neumann boundary conditions 



u x (0) = u x 



= 0, v x (0) = v x 



(3.10) 



the system (|2.20J) can be written in the form (|3.5|) where 
C\ is symmetric and Cq 1 symmetric and positive definite. 
As before, the minimum eigenvalue of (13.51) is given by 
the minimum of the Rayleigh quotient ()3.6j) where this 
time the minimum is evaluated over all u(X) satisfying 
the Neumann conditions p.lOj) . 

The operator C\ has a negative eigenvalue 



Mi 



with an eigenfunction 



Vfc 4 + 1 - k 2 
2~k 2 



dn 2 



X 



V2^F 



zi{X) 



which satisfies (|3.9[) . Therefore, the minimum of the 
Rayleigh quotient is negative on the space of u satisfy- 
ing (|3.9p ; the vector eigenvalue problem (|2.20p has a real 
eigenvalue A, and the growth rate (I2.22p is positive. 

Recalling that the eigenfunction y(X) can be chosen 
either even or odd, the condition 2/x(0) = selects the 
even function. The boundary conditions (13.101) imply 
then that yx(—L/2) = and so the eigenfunction y(X) 
is L-periodic: 



V \-2 =V 2 



L 



L 



Vx\- 1 =VX\ g 



L 



(3.11) 



We have thus established instability of the dn~ solu- 
tion against perturbations periodic with the period L^ n , 
the period of the dn cnoidal wave. The dn~ wave is un- 
stable for any choice of the parameters h and 7, and any 
k. 



C. Instability of dn + : antiperiodic perturbations 

So far we have demonstrated the instability of the cn~ 
and dn - waves, with the unstable perturbations exhibit- 
ing a monotonic growth. Another solution that turns out 
to be prone to the instability of a similar type is the dn~ 
wave; however this time our proof will only be valid in a 
part of the (h,j,k) parameter space. 

In the case of the dn + wave, the operator Co has an 
eigenvalue Ai = 1/(2 — k 2 ) — £ with an eigenfunction 



yi(X) = sn 



X 



.V2^W 

The eigenfunction y\ satisfies mixed boundary conditions 

(3.12) 



2/(0) = 0, y x ( - ) = 



and does not have zeros inside the interval (0,L/2). 
Hence the eigenvalue Ai is the smallest eigenvalue of Cq 
under the boundary conditions (|3.12l) . 

The eigenvalue Ai is positive and the operator Cq is 
positive definite in two adjacent parameter regions. One 
parameter region is 



1 



(3.13) 

with k taking any values between and 1. The second 



region is 



y I + 7 2 < h < v/TT 



r 



(3.14a) 



with the elliptic modulus being bounded from below: 



k z > 



3 1 1 



2 2 ^h 2 - 7 2 



(3.14b) 



We now assume that the parameter vector (h, 7, k) lies 
in one of the above two regions. 

The operator C\ has a negative eigenvalue 2(fc 2 — 
l)/(2 — fc 2 ) with an eigenfunction satisfying the boundary 
conditions (|3.12|) : 



z\{X) = sn ( — , k ) dn [ — , k 



V2~W 



V2^W 



Therefore the minimum of the Rayleigh quotient (|3 .€>[) on 
the space of functions satisfying (|3.12l) is negative and the 
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symplectic eigenvalue problem (I2.20|) has a real eigen- 
value A. This means that the dn + wave is unstable to 
perturbations satisfying 



u (0) = u x 



= 0. -'(Hi = r A ( | I = U. 



The boundary condition y(0) = singles out the odd 
eigenfunction; hence it satisfies the antiperiodicity con- 
ditions on the interval (— L/2,L/2): 



y 



-y 



yx 



-yx 



We conclude that when h, 7 and k belong to the region 
(I3.13[) + (|3.14p . the cnoidal wave dn + is unstable under 
perturbations of period twice the period of the wave. 



IV. THE UNPERTURBED NONLINEAR 
SCHRODINGER: STABILITY OF THE cn WAVE 



When h = 7, Eq. (|2.10p gives £ — and the spectrum 
of symplectic eigenvalues (12.20)) coincides with the spec- 
trum of stability eigenvalues of the cnoidal wave of the 
undamped undriven NLS. 



A. Small eigenvalues: general setting 

In this subsection we will obtain small symplectic 
eigenvalues. 

Factorising the eigenfunction y into a periodic func- 
tion y and an exponential, y — ye lKX , the symplectic 
eigenvalue problem ()2.20|) becomes 



ny = 2i K y x - n 2 y + \jy. 



(4.1) 



Without loss of generality, the period of the function y 
could be taken equal to the period of the potential of the 
operator %: L cn in the case of the cn- wave, and L& n in 
the case of the dn-w&ve. This will indeed be our choice 
in the dn situation. However, in the case of the cn wave, 
it is convenient to regard y as a 2L cn -periodic function 
- that is, choose the period of y to coincide with the 
period of the cnoidal wave. This convention is equivalent 
to the previous one, and is equally general. 

The symplectic spectrum includes a four-fold zero 
eigenvalue; associated with these are two periodic eigen- 
vectors and two generalised eigenvectors. When k is 
small, we expand the eigenvalues and eigenfunctions in 
powers of n: 

y = % + k& + n 2 y 2 + 

X = X 1 k + X 2 k 2 + ... (4.2) 

Substituting in (|4.1I) we equate coefficients of like powers 
of K. 



The coefficient of k° gives Wyo = 0. The general solu- 
tion is of the form 



^0 = 



( M 



x 



(4.3) 



where A and B are two arbitrary constants. Next, the 
order k 1 produces 

uy x = 2id x % + (4.4) 

or, componentwise, 

Cim = 2iAd 2 x q - XiBq, (4.5) 



£ vi = {2tB + X 1 A)q x , 



(4.6) 



where u\ and v\ are the top and bottom components of 
the vector y%: $\ = (ui,vi) T . We note that 

dd x (Xq) = -2q, 
d{Xq x ) = -2d\q, 
£ (Xq) = ~2q x ; 

therefore one solution to Eq. ([4.5[) is 

ux = ~iAXq x + ^-B(Xq) x , (4.7) 
and one solution to (14.61) is 



(4.8) 



vi = --(X 1 A + 2iB)Xq. 
However neither of these solutions is periodic. 



B. cn wave 

From this point on, our analysis depends on which peri- 
odic solution we consider. We start with the cn wave. In 
this case both components of the zero-order approxima- 
tion y are 2L-periodic, and we will attempt to construct 
its perturbation y with the same period. Here L is our 
short-hand notation for L cn . 

In order to obtain a periodic ui, we can add to u\ a 
multiple of d^q, the nonperiodic homogeneous solution of 
equation (|4.5I) : 



dkq(x) = sn(£, fc)dn(£, k) 
E(a,m£,k) X 



(2k 2 - 1)V2 



(l-fc 2 )(2fc 2 - 1) (2fc 2 -l) 2 

1 k 2 

' sn 2 (C,fc) 



2fc 2 - 1 1 - k 2 



(4.9) 



In l[49]) . £, = X/V2k 2 - 1 and E(ip, k) is the incomplete 
elliptic integral of the second kind: 

E{ip,k)= [ V \fl - k 2 sin 2 OdO. 
Jo 
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FIG. 1. Spectra of the cn wave with h = 7. Top row: the spectrum of the symplectic operator J~ 1 H with £ — 0. This is 
essentially the spectrum of the cn-wave solution of the unperturbed NLS. The dashed line in (a) and (b) makes the angle <j> 
to the horizontal axis, where <j> is as in (|4.21|) . Bottom row: the corresponding spectrum of the linearised operator (|2.18|l with 
h = 7. In these plots, 7 = 0.5. 



The resulting nonhomogeneous solution 



Ui = -lAXqx + Y B ( X( l)x + Cd k q (4.10) 

is even and so the periodicity condition u\(L) = ui(—L) 
is satisfied automatically, for any C; therefore, one only 
needs to satisfy dxUi(L) = dxUi(—L). For the odd 
dxUi(X), this reduces to 

dxu^L) = 0. (4.11) 

Substituting (j4~T0f in (jlTi) we find 



C = (\ 1 B-iA)—, 
J->k 

where L k = d k L. 

I 



At the order k 2 we obtain 

Ciit2 = 2idxu\ — u — X2V0 — X1V1, (4-14) 
C0V2 = 2idxvi — vq + X2U0 + Aiiti, (4-15) 



In a similar way, we take 

Vi = -(^\ 1 A + iBjXq + Dz, (4.12) 

where z is the nonperiodic homogeneous solution of (14.61) , 

z(X)=£cn(£,fc) 

Eq. (|4.12[) is odd and its derivative dxv\ is even; hence 
the periodicity condition dxVi(L) = dxV\{—L) is in 
place and we only need to make sure that v\{L) — 
V\{—L). This amounts to 

vi(L)=0. (4.13) 

Substituting Eg. (14.121) in (|4.13[) . the constant D is iden- 
tified: 

D = -(\ 1 A + iB)— . 



I 

where (u2,V2) T — 3^2- The solvability conditions for 
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these two equations are, respectively: 

2i(q x \d x u 1 ) - A(q x \qx) - Mlx\vi) = 0, (4.16) 
2i(q\d x vi) - B(q\q) + Xi{q\ux) = 0. (4.17) 



Here 



(f\9) = f f{X)g{X)dX. 



Substituting for u\ and v\ we obtain a system of linear 
equations 



M 



A 



mnX 2 + fin imi 2 \i 



A 



0. 



B J y im 2 i\i m 22 \\+n 22 J \ B 

(4.18) 

where the elements of the matrix M. are defined by 

™n = iMxWi, nn = (qx\2px - qx), 
mu = (qx\qx +Px +w), m 2 i = -(q\p + w x ), 



m 2 2 = -{q\q +p), n 22 = (q\2w x - q). 



(4.19) 



Here 



L K 
p = Xq x + —qk, w = Xq+—z. 



The integrals are easily evaluated: 



mn = 



fc(l-fc 2 ) (LN) k 
2(2k 2 -l) N 

mi2 = m 2 i 



nn 



2k L 2 
(2fc 2 -l) 2 L? 
E(E - K) 



{2k 2 - 1)(1 - fe 2 ) K k L k ' 



m 22 = - 



1 (LN) k 



2 L k 



-, n 22 



4k K 2 
(2k 2 



where 



N = J q 2 (X)dX = 



4(1 -k 2 ) 
V2k 2 - I 



and i given by Eq. (|2.6[) . Setting det.M = yields a 
biquadratic equation 



where 



and 



aoAi + b a \j + c = 0, 



[frml 2 , 8kL 2 N 2 L 



(4.20) 



16k 2 L 4 



Co 



(2fc 2 - l) 2 ' 

The roots of the equation are given by ±Ai,±A*, where 



Al=P ,«, „ t , = ^n^_£. (4 . 21) 



The roots have nonzero real parts; this implies that 
the cn-solution of the unperturned nonlinear Schrodinger 
equation is unstable for all k. This fact is known in lit- 
erature. 

The eigenvalues A = AiK + 0(k 2 ), with small k, lie on 
four rays emanating out of the origin on the complex-A 
plane. The full spectrum of symplectic eigenvalues, ob- 
tained numerically, is displayed in the top row of FigfT] 
The result (|4.21[) for the dominant behaviour of the eigen- 
values near the origin is shown by the dashed line. The 
analytical result is seen to accurately reproduce the nu- 
merically computed eigenvalues. The four rays bend and 
join, pairwise, forming an eight-shaped curve centred at 
the origin (FigfT] (a-c)). In addition, the spectrum fills 
the imaginary axis of A (with a gap). 



C. The unperturbed nonlinear Schrodinger: 
Stability of the dn wave 

In the case of the dn wave, the eigenvalue problem 
(|4.1[) and expansion (|4.2[) remain in place. The null eigen- 
function J^o is periodic with the period of the dn wave: 
yo(X+L) = y (X), where L = L<jn is as in ()2.7[) . and we 
will assume that the eigenfunction y satisfies the same 
boundary conditions. 

The e^corrections ui and vi are still given by equa- 
tions (|4.10p and (|4.12j) . respectively, where the homoge- 
neous solutions are, this time, 



dkq 



and 
where 



fcsn(£, fc)cn(£, k) 
(2-fc 2 )V2 



£(am() X 
1-k 2 ~ (2 - k 2 ) 

[(2-fc 2 )cn 2 (£,fc)-l] 



(l-fc 2 )(2-fc 2 ) 
fc 2 cn(£, fc)sn(f , k) - 2dn(£, fc)£(am£, k) 
X 



V2~W 



Since the function 



Ai 



Ui(X) = -lAXqx + -^B{Xq) x + Cq, 



(4.22) 



is even, the only periodicity condition that needs to be 
verified, is 



0\ "j ( - 



0. 



(4.23) 



Substituting (j4.22|) in (|4.23[) , the constant C is evaluated 
to be 



C = (AiB - iA) 



L 
Lk 
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where L k — d k L. On the other hand, the function 



Substituting for q we get 



Vi = -[ -AXi +iB ) Xq + Dz 



is odd; hence the periodicity condition reduces to 



This gives 



D = -(\ 1 A + iB)-. 



At the order k 2 we obtain equations (|4.14[) - (|4.15l) . with 
the solvability conditions (I4.16p ~ d4.17p . where, this time, 
the scalar product is defined as 



L/2 

(/| fl > = / f(X)g(X)dX. 

-L/2 



The elements of the matrix (14.181) are given by Eqs. (|4.19p 
where, this time, the functions p(X) and w(X) are given 
by 



mn 



k(l - fc 2 ) (LN) k 
4(2 -fc 2 ) 2 N 



-, "-11 



mi2 = m 2 i = - 

x (yi - fc 2 iv) 

1 (LN) k 

m 2 2 = 7: — z , n 2 2 = 

l L k 

where L is given by (|2.7I) and 



fc 2 L 2 
= {2-k 2 ) 2 L k 
k' 2 L 
(2 - fc 2 ) 2 NL k 
L 



1-fc 2 L 2 
_ (2 — fc 2 ) 3 TV' 



/•L/2 2 

A = / q 2 (X)dX = _ F. 

J L/2 



(4.24) 



Hence we find the coefficients of the biquadratic equa- 
tion (|4~20l) : 



a a = [(LN)k\ , c 



8k 2 r4 



(2 - fc 2 ) 3 



and 



bo = - T77. To^ LN 2 



k(2-k 2 ) \Nj k 
The discriminant of the equation is 



V 



k(2-k 2 ) 



T (k' 2 L k N - N k L)(L k N - k' 2 N k L). 



p = 2Xqx + -^q k (X), w = Xq+§z(X). 
L k hj 



Since L k > and N k < 0, this is positive and so Eq. (|4.20p 
has two positive roots, (A 2 ) > and (A 2 )& > 0. 



The conclusion is that the dn-wave solution of the 
unperturbed NLS is unstable, for any k. This fact is 
known to workers in the field. Furthermore, the real pos- 
itive symplectic eigenvalues A translate into real positive 
growth rates 77. This implies that the dn wave is also 
unstable as a solution of the damped-driven NLS with 
h = 7. 



V. THE DAMPED-DRIVEN cn WAVES WITH 

h = 7 

That the cn wave of the unperturbed NLS is unsta- 
ble for any k does not mean, however, that the damped- 
driven wave with h = 7 is necessarily unstable. In the 
damped-driven situation, the stability is determined by 
the linearised eigenvalues rj, not symplectic eigenvalues 
A. Transforming A to ?y by the rule (|2 . 19[) , the symplectic 
spectrum shown in the top row of FigQ]is mapped to the 
spectrum shown in the bottom row of the same figure. 



A. Stability to long-wavelength perturbations 

The symplectic spectrum includes a double zero eigen- 
value resulting from the phase invariance of the unper- 
turbed NLS equation (invariance w.r.t. constant and 
constant-velocity phase rotations) and another double 
zero resulting from its translation and Galilian imparl- 
ances. The map (|2.19l) leaves only two eigenvalues at the 
origin: one resulting from the translation symmetry and 
the other one corresponding to the symmetry w.r.t. the 
velocity boosts. 

Next, for small k <C 7, the small-77 branch of the map is 
simply rj = (27) _1 A 2 . Hence any ray A = e*^r emanating 
out of the origin on the A-plane and making the angle <p 
to the real axis, is mapped to a ray making double that 
angle to the real axis on the 77-plane: 77 = e 2l ^r'. The 
branches of the spectral curve with Re A > emanate 
out of the origin as two rays making the angles ±0 to 
horizontal. Eq. (14.211) implies that the value of cot cf) is 
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FIG. 2. Spectra of the dn wave with h = 7. Top row: the spectrum of the reduced operator (|2.20[) with £ — 0. Bottom row: 
the corresponding spectrum of the linearised operator (|2.18[1 with h = 7. In these plots, 7 = 0.5. 



smaller than 1 for all k > -j= (see the Appendix) and 
so cf> lies between j and ^. These two rays (and their 
negatives) are mapped to rays with Re 77 < 0. Therefore, 
no instability of the damped-driven cnoidal waves with 
h = 7 is associated with the neighbourhood of r\ — 0, no 
matter what is the value of k. 



This conclusion is valid for all 7, including very small 
ones. A natural question, therefore, is how the stable 
77-spectrum becomes unstable when 7 reaches zero; that 
is, how can a diagram from the bottom row in Fig IT] (two 
rays at the angle 2<fi to horizontal) evolve into the corre- 
sponding diagram in the top row (two rays at the angle 
4> to horizontal). The answer is that for small 7, 7 « k, 
the map (|2.19p reduces to 77 = A. On a relatively large 
scale, namely on the scale |?7| ^> 7, the spectrum of r\ 
looks indistinguishable from the corresponding spectrum 
of A; it includes an eight-shaped curve centred at the ori- 
gin. However if we zoom in on very small 77 (\r}\ ~ 7), 
we will observe that the four rays do not reach the origin 
on the 77-plane, but form a shape shown in the inset to 
FigEl As 7 — > 0, the box shown in the inset shrinks to 
the origin. 



B. The cn waves with h = 7: arbitrary 
perturbations 

The periodic wave is unstable if the linearised- 
spectrum curve crosses into the Re r\ > half-plane. For 
the crossing points Re r\ = 0, Eq. (|2.19l) gives 



ImA = ± 



Re A 



VI -(Re A/7) 2 ' 



(5.1) 



This is an equation for two curves on the A-plane. These 
two curves are straight lines at the origin, making 45° to 
the horizontal. As ImA — > ±00, the curves are asymp- 
totic to the vertical straight lines Re A = ±7 (see FigH]) . 
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FIG. 3. Spectrum of linearised eigenvalues of the cn solution 
with very small 7. (In this plot, 7 = 0.005 and k 2 = 0.8.) 
On a large scale, the spectrum of 77 looks like the spectrum 
of A; however zooming in on the neighbourhood of the origin 
(dotted box, enlarged in the inset) the shape typical for 77- 
curves emerges. 




FIG. 4. The symplectic spectrum of the cn wave (solid) and 
the curves (|5.1|l (dashed). For small 7, the curves intersect the 
spectrum. The "internal" section of the spectrum corresponds 
to linearised eigenvalues with Re 77 < whereas the part which 
lies outside the dashed curves is "unstable": Re 77 > 0. For 
large 7, the entire symplectic spectrum is stable. In this plot, 
k 2 = 0.65. 



For the given h = 7, the cnoidal wave is unstable if the 
corresponding pair of curves (|5.ip intersects the locus of 
the symplectic spectrum on the A-plane. Since the rays 
of the eight-shaped spectral curve make the angle greater 
than 45° to horizontal at the origin, the curves (|5.1I) do 
not have to intersect the spectral curve. They will cross 
through the "eight" only when 7 is large enough (FigQ} . 
For some critical value 7 = 7 c (fc), the curves will just 
touch the "eight". This critical value defines the range 




FIG. 5. The critical value of the elliptic modulus as a function 
of 7. Shaded is the stability domain; for k < k c (-f) the wave 
is unstable. 



of stability of the cn wave for the given k: the wave is 
stable when 7 > j c (k) and unstable otherwise. 

Note that when h — 7, the original and scaled damping 
coefficients coincide: 7 = 7. Therefore, on the (7, h)- 
plane, the stability region is given by the ray h = 7 with 
7 > 7 C (k). The inverse function k c (-j) [where k c (j c (k)) = 
k] defines the boundary of the stability region on the 
(7, fc)-plane (FigJS]). 



VI. SOLUTIONS cn+ AND crT : h > 7 

In this section, we calculate the deformation of the 
linear spectrum associated with the cnoidal waves cn + 
and dn + as h grows from the value h = 7. 

It is not difficult to evaluate the change of the symplec- 
tic spectrum as h deviates from 7 and hence the param- 
eter £ deviates from zero. The corresponding spectrum 
coincides with the spectrum of the undamped NLS with 
nonzero parametric driving. 

In addition to expanding y and A as in (|4.2[) , we let 
£ = £qk 2 in the operator (|2.9a[) [which, in turn, is a ma- 
trix element of the operator (|2.2ip j. Substituting these 
expansions in equation (|4.1|l . the orders k° and k 1 pro- 
duce the same equations as in section ITVl with the same 
set of solutions. The first difference from the case £ = 
arises at the order k 2 where the equation (|4.15|) is re- 
placed with 

C v 2 = 2id x vi - (1=f£o)«o + A 2 uo + Aiui, 

and the corresponding solvability condition (|4.17p with 

2t(q\d x v 1 ) - B(l T £o)(q\q) + Ai(<?K) = 0. 



[On the other hand, the solvability condition (|4.16l) is 
not changed.] As a result, the only matrix element of the 
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matrix A4 that is altered by allowing a nonzero £, is ni2 
which becomes 



Ak K 2 



TL22 



±£ N. 



y/2k 2 - 1 K k 
The quartic equation (|4.20p is replaced with 

aoA? + (b T ^q6i)Ai + (co T £qCi) = 0, 



where 



6i = -2(LN) k L k 



ALN 2 



ci = 



(l-fc 2 )(2fc 2 -l) 3 



are two negative constants. 

Defining such that Ai = 
becomes 



1/2 

£ /Lt, the quartic equation 



6 ° n-l, 



(6.1) 



Here £q is a positive parameter. 

Consider the top sign in (|6.ip . that is, consider eigen- 
values pertaining to the cnoidal wave cn + . For £ a = oo, 
two roots of equation (|6.1[) are equal to zero, fi = 0, while 
the other two roots take opposite pure imaginary values: 
/i = ±i(— 61/ao) 1 ^ 2 . As £q decreases from large values, 
two pairs of opposite roots move along the imaginary 
axis, collide pairwise, and appear into the complex plane, 
forming a quadruplet ±/x, ±/x*. As £q — > 0, the four com- 
plex roots diverge from the origin along the straight lines 
making the angles ±0 to the horizontal, with <fi as in 
(I4.21j) . The trajectories of the roots of the equation (16. ip 
as £0 is varied are shown in Fig|6l 

The symplectic eigenvalues A are related to the roots 
\i by the scaling A = S 1 ^ 2 ^. Thus for any small £ , the 
small- A part of the symplectic spectrum of the cn + wave 
is obtained by scaling down Fig|6l 



Appendix A: Some useful inequalities involving 
elliptic integrals 

1. First we show that the number of particles captured 
in one period of the dn wave, eq. (|4.24[) . is a monotonically 
decreasing function of k. A quick way to establish this is 
to note that the derivative 



N k = 



[2E - (2 - k 2 )K] 



fc(2-fc 2 ) 3 / 2 ' 
can be related to the manifestly positive integral 

/' •^K' t) *->-"-W 

Jo dn z (£,k) k 2 

2. Next, we prove that the derivative 

Tk^ = W^) [E2 - (1 - k2)K2] (A1) 
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FIG. 6. Solid curves: trajectories of the four roots of Eq. (|6.1[) 
as £0 is varied from to 00. The dashed lines are plotted at the 
angles ±cj> to the horizontal, with <f> as in (|4.2ip . Dots: numer- 
ically obtained symplectic eigenvalues. In this plot, k 2 = 0.8. 
£ = 0.01 



with N as in (|4.24|) , is positive. To this end, it is sufficient 
to observe (a) that the expression in the right-hand side 
of (|A1[) equals zero at k = 0, and (b) that the derivative 
of this expression is positive: 



dk 



[E 2 — (1 — k 2 )K 2 ] = -(E-K) 
k 



3. Since (LN) k > while N k < 0, one is led to con- 
clude that L k > 0, the period of the dn wave is a mono- 
tonically growing function of k. An independent way to 
see this is to note that the derivative 



dL dn 2V2^F 



E 



2K 



dk k \l-k 2 2-k 2 / 

can be related to the manifestly positive integral 

f K 2^,i 2„,n„ (2 - k 2 )E - 2(1 - k 2 )K 
J sn 2 (£, fc)cn 2 (£, k)d£, = y - >-^± >—. 

4. Finally, we prove the inequality 
k' K(k) - E(k) 



k 



E(k) 



< 1 



(A2) 



for all i < k 2 < 1. Here k' — yl — k 2 is the complemen- 
tary modulus of the elliptic integrals. 
We start by writing (|A2[) as 



g(k) < 1, 



(A3) 



where we have defined 
9(k) 



k 1 dE 

E{k)~d~k' 
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At the point fco — 1/v2j the integral K(k) has the 
property dK/dk — K which allows to relate K and E, 



K 



fc' 2 (l + fc) 



-E, 



and subsequently evaluate g(k ): 

a/2-1 



g(ko) 



\/2 + l' 



Therefore, at the point fco the value of the function g(k) 
is smaller than 1 and will remain smaller than 1 in some 
neighbourhood of this point. We will now show that this 
neighbourhood extends all the way to fc = 1. 

Assume the contrary, that is, assume that there is a 
point fc», fco < fc* < 1, such that 

g(K) = 1. (A4) 

Using the hypergeometric equation satisfied by E(k), 



we obtain a Riccati equation 

df_ = P 
dk k 

for the function 



fc' 2 



(A5) 



/(fc) = -g(k). 
Integrating both sides of (| A5|) from fco to fc* we obtain 



/(*♦)- /(*o) 





7 2 " 


/fco 


. k + fc' 2 . 



dfc. 



(A6) 



Since /(fc) < fc/fc' in the interval (fco,fc*), the integral 
admits a simple bound: 
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k 



k 

fc 72 



dfc < 



1 / fc 
fc I fc 7 



_fc_ 

fc 72 



dfc. (A7) 



The integral in the right-hand side of (|A7|) is tabular and 
hence (|A6j) gives 



/(fc*) < 1 + ln 



fc 2 



i-fcr 



(A8) 



where we have also used /(fco) < 1. 

On the other hand, one can readily show that the func- 
tion 



Fi(fc) = 1 + ln 



fc 2 
o 



1 



is smaller than F 2 {k) = k/k' for all fc > fco- [Indeed, we 
have Fi(fc ) = F 2 (fc ) = 1 but dFi/dk < dF 2 /dk for all 
< fc < 1.] In particular, Fi(fc*) is smaller than ^(fc*), 
that is, 



1 + ln ^ 



< 



1 - fc 2 fc' ' 



(A9) 



where k'„ = yjl - fc 2 . Using (|A"9|t . the inequality (fMJ 
becomes /(fc*) < fc*/fc* or, equivalently, 



(A10) 



which contradicts (|A4|) . 

The contradiction proves that no point fc* < 1 satisfy- 
ing (|A4I) can exist. Therefore the inequality (|A3|) (and 
thus, (|A2[) ) remain valid for all fc between fco and 1. 
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